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ABSTRACT 

Context. Protoplanetary nebulae typically present non-spherical envelopes. The origin of such geometry is still controversial. 
There are indications that it may be carried over from an earlier phase of stellar evolution, such as the AGB phase. But how 
early in the star’s evolution does the non-spherical envelope appear? 

Aims. Li-rich giants show dusty circumstellar envelopes that can help answer that question. We study a sample of fourteen 
Li-rich giants using optical polarimetry in order to detect non-spherical envelopes around them. 

Methods. We used the IAGPOL imaging polarimeter to obtain optical linear polarization measurements in V band. Foreground 
polarization was estimated using the field stars in each CCD frame. 

Results. After foreground polarization was removed, seven objects presented low intrinsic polarization (0.19 — 0.34)% and two 
(V859 Aql and GCSS 557) showed high intrinsic polarization values (0.87 — 1.16)%. This intrinsic polarization suggests that 
Li-rich giants present a non-spherical distribution of circumstellar dust. The intrinsic polarization level is probably related to 
the viewing angle of the envelope, with higher levels indicating objects viewed closer to edge-on. The correlation of the observed 
polarization with optical color excess gives additional support to the circumstellar origin of the intrinsic polarization in Li-rich 
giants. The intrinsic polarization correlates even better with the IRAS 25 pm far infrared emission. Analysis of spectral energy 
distributions for the sample show dust temperatures for the envelopes tend to be between 190 and 260 K. We suggest that dust 
scattering is indeed responsible for the optical intrinsic polarization in Li-rich giants. 

Conclusions. Our findings indicate that non-spherical envelopes may appear as early as the red giant phase of stellar evolution. 

Key words, polarization - stars: circumstellar matter - infrared - stars: atmospheric parameters - stars: abundances - stars: 
late type giants - stars: evolution - element: Li 


1. Introduction. 

Mass loss plays a central role in the late stages of stellar 
evolution. In particular, asymptotic giant branch (AGB) 
stars may be responsible for up to 60% of the (all stars) 
interstellar dust input into the interstellar medium (Gehrz 

H2EB9- 

Our current understanding of stellar evolution places 
the AGB stars as precursors of protoplanetary nebulae 
(PPN). On the other hand, we know from direct imag¬ 
ing that PPN typically present non-spherically symmetric 
envelopes. While extrinsic (to the PPN) effects such as 
binarity have been said to cause the non-sphericities (for 
a review see Balick & Frank 2()(')2 i. could it not be instead 
that the asymmetry is carried over from the earlier, AGB 
phase of stellar evolution? 


Send offprint requests to: A. Pereyra, e-mail: 
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* Based on observations obtained at the Observatorio do Pico 
dos Dias , LNA/MCT, Itajuba, Brazil. 


One piece of evidence showing that departures from 
spherical symmetry exist around AGB stars is that light 
from these objects may show some degree of linear polar¬ 
ization (Coyne & Magalhaes 119771 ft 0701 Magalhaes et al. 
I1986alll986bl Kahane et al. 119971 Magalhaes & Nordsieck 
mm- Further evidence of asphericity in the envelopes 
of late-type giant and supergiant star envelopes comes 
from the details of OH maser emission profiles (Collison 
& Fix 11992 1). KI resonant scattering (Plez & Lambert 
11994 ?), rings of SiO maser emission (Diamond et al. 119941 
Greenhil et al. and OH radio images (Chapman et 

al. rHMll . 

Aspherical symmetries, such as those found by 
Trammell et al. (1994) in post-AGB stars from spec- 
tropolarimetry, may then be understood naturally, since 
such symmetries are already present in earlier evolution¬ 
ary stages. The observed and more obvious non-spherical 
symmetries in PPN (Balick & Franck 120021) and the in¬ 
ferred assymetries in other evolved objects obtained from 
polarimetry (Johnson & .)one> : 1991 Parthasarathy & Jain 
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1 992) Parthasarathy et al. 2DTF5 ) are also consistent with 
the origin of the asymmetries early in the AGB phase. 

Aside from the fact that polarimetry per se allows the 
study of otherwise unresolved objects, the question arises 
as to how early in stellar evolution the non-spherical sym¬ 
metry appears. In this paper we use the fact that Li-rich g 
red giants (RG) show dusty circumstellar envelopes to ex- § 
plore their environments and to help answer that question. 

Over the past two decades around 40 red giants were 
found to have Li abundances that are 100 times larger 
than the mean values observed in red giants (Brown et al. 
Mi de la Reza & Drake H2S3- This number indicates 
that about 2% of the Population I red giants show signif¬ 
icantly larger lithium abundances than expected by dilu¬ 
tion due to mixing by classical convection. In some of the 
giants, the lithium abundance reaches values that are sim¬ 
ilar to (and even larger than) the Pop I value (meteoritic, 
open clusters, etc.), around logN(Li) = 3.3. Castilho et al. 
( 1999a ) obtained very low Be abundances for two Li-rich 
giants (LRG) providing evidence that the original Li in 
these stars must have been almost completely destroyed 
and that the high Li abundances in the Li-rich red giants 
are due to Li production in these stars. 

It seems to be clear that LRG are quite normal stars, 
except for their high Li abundance and large infrared ex¬ 
cess (Castilho et al.jTSSSI de la Reza et al. lUWbt Castilho et 
al. roi) . The similarities to normal giants, namely mass, 
chemical composition, temperature, and metallicity, com¬ 
bined with the far-infrared emission, indicate that LRG 
do not form a unique class of objects but are ordinary 
low-mass stars observed during a short phase of their evo¬ 
lution, when Li is created. 

If indeed all low mass red giants go through a phase of 
Li production in the RGB (that could be cyclic), together 
with an increase in mass loss seen in the IRAS color dia¬ 
gram (Gregorio-Hetem et al. 119931 de la Reza et al. 1199611 . 
and since some of them reach a larger lithium abundance 
than the Pop I abundance, they could be an important 
source of Li enrichment in the Galaxy. Some information 
needed to quantify the LRG contribution for the interstel¬ 
lar medium Li enrichment remains unknown, such as: the 
maximum Li abundance reached for each star and its re¬ 
lation with stellar mass and/or metallicity, the duration 
of the Li production phase(s), the mass loss process, and 
rate, and the simultaneity of the envelope ejection with 
the Li production. 

Up to now only one envelope of LRG has been stud¬ 
ied in detail. The ~ 3x4 arcmin envelope of HD 65750 
(Castilho et al. 1 99 m has a butterfly geometry and the 
present mass loss rate of the star is not enough to form the 
observed envelope. Witt & Rogers mm proposed that a 
past and more efficient mass-loss event about 32,000 years 
ago was responsible for the observed structure. 

The aim of this work is to measure the polarization 
of LRG stars so as to study the spatial distribution of 
the circumstellar dust around these objects. Scattering of 
the stellar radiation by the circumstellar dust can pro¬ 
duce polarization that is measurable in objects with non- 
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Fig. 1 . Histograms of (a) the observed and (b) intrinsic 
polarization of LRG. The NG in the sample are indicated 
in black. 


resolved envelopes. The amount of polarization will de¬ 
pend upon the density and nature of the scatterers, as well 
as on the aspect angle and exact geometry of the envelope 
nvia.galhaes tl 992! If intrinsic polarization is detected, we 
can infer that an asymmetric spatial distribution of cir¬ 
cumstellar dust is present in these objects. 

Here we present the results of optical linear polariza¬ 
tion measurements of fourteen LRG and two normal giants 
(NG). In Sect. 2 we describe the observations and data re¬ 
duction along with the calculations of foreground and in¬ 
trinsic polarizations. In Sect. 3 we show the stellar param¬ 
eters of our sample and a general discussion is presented 
including correlations between the polarimetric data and 
optical and near infrared excess color. Correlation with 
IRAS colors also are explored and spectral energy distri¬ 
butions used to estimate the dust temperature associated 
to the circumstellar envelopes. The conclusions are drawn 
in Sect. 4. 

2. Observations and data reduction 

The observations were made using IAGPOL, the IAG 
imaging polarimeter (Magalhaes et al. 119961) . at the f/13.5 
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Table 1. Log of observations. 


Object 

Log e(Li)“ 

l 

b 

V" 

Wav. pos. 

IT 

Aper. 

Date 

ID 


(dex) 



(mag) 


(sec.) 

(") 



(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

Li-rich giants 

HD 90082 

0.20 (0.15) 

286.27 

-3.91 

7.50 

16 

8 

5.5 

06/20/98 

a 

HD 95799 

3.10 (0.15) 

289.19 

1.17 

8.14 

16 

12 

5.5 

06/20/98 

b 

HD 96195 

0.40 (0.20) 

291.12 

-2.57 

7.94 

16 

13 

6.7 

06/20/98 

c 

HD 120602 

1.90 (0.15) 

338.54 

64.21 

6.01 

16 

1 

3.7 

06/20/98 

d 

PDS 68 

2.00 (0.40) 

315.72 

19.10 

12.80 

4 

300 

4.3 

06/20/98 

e 

HD 146850 

1.60 (0.15) 

359.59 

24.44 

6.10 

6 C 

1 

3.7 

06/20/98 

f 

GCSS 557 (V385 Set) 

0.70 (0.15) 

17.10 

-1.37 

13.26 

4 

300 

3.7 

06/20/98 

g 

HD 176588 

1.60 (0.15) 

30.08 

-4.18 

6.89 

16 

2 

4.9 

06/20/98 

h 

IRAS 19012-0747 

2.60 (0.15) 

27.46 

-6.28 

11.17 

16 

60 

3.1 

06/20/98 

i 

IRAS 19038-0026 

0.60 (0.25) 

34.32 

-3.49 

14.48 

8 

300 

3.7 

06/20/98 

j 

HD 178168 

0.90 (0.15) 

32.54 

-4.72 

9.06 

16 

20 

5.5 

06/20/98 

k 

HD 112127 

2.80 (0.15) 

0.96 

89.34 

6.87 

16 

1 

4.9 

06/21/98 

i 

V859 Aql (PDS 100) 

2.50 (0.15) 

42.29 

-6.28 

10.44 

16 

40 

4.9 

06/21/98 

m 

HD 203251 

1.40 (0.15) 

35.57 

-39.97 

8.00 

16 

8 

7.3 

06/21/98 

n 

normal giants 

HD 190664 

<0.00 (0.30) 

37.92 

-18.48 

6.47 

16 

1 

3.7 

06/21/98 

o 

HD 124649 

<0.00 (0.30) 

315.53 

7.45 

7.86 

16 

6 

6.1 

06/21/98 

p 


Notes: (a) log e(Li) were obtained from Castilho et al. nr except for: HD 120602 and HD 112127 (from Brown et al. 11989!) . 

HD 146850 (from Castilho et al. 11995 ). V959 Aql (from Reddy et al. 12002 3. and HD203251 (from Fekel & Balachandran). 
The value for PDS 68 is an estimation based on the equivalent width of the Li I (6707.8A) line (from Gregorio-Hetem et al. 
119921 . and the values for HD 190664 and HD 124649 are higher limits (Castilho et al. T999ai l. The errors are indicated in 
parenthesis; (b) V were obtained from Castilho ( 1999bt . except for: HD 120602, HD 112127, HD 203251, and HD 190664 (from 
GCPC, Mermilliod et al. 11997 1. PDS 68 and V859 Aql (from Gregorio-Hetem et al. 11992^ . HD 124649 from CDS, and IRAS 
19038-0026 (from the GSC 2.2 catalogue, using magnitude in the visual phographic band as V); (c) 8 waveplate positions were 
initially observed but two of them were saturated. 



Fig. 2. Correlation between observed polarization and 
galactic latitude. The black dots are the LRG, and the 
white dots are the NG. The letters indicate the ID for 
each object as it appears in column 10 of Tabled 

Cassegrain focus of the 0.6m IAGUSP Telescope at the 
Observatorio do Pico dos Dias (OPD), operated by the 
Laboratorio Nacional de Astrofisica (LNA), Brazil. When 
in linear polarization mode, the polarimeter incorporates a 
rotatable, achromatic half-wave retarder followed by a cal- 
cite Savart plate. This provides two images of each object 
in the field, separated by 1mm (corresponding to 25.5" at 


the telescope’s focal plane), with orthogonal polarizations. 
One polarization modulation cycle is covered for every 90° 
rotation of the waveplate. The simultaneous imaging of 
the two beams allows observation under non photomet¬ 
ric conditions, and at the same time the sky polarization 
is practically canceled. The polarimeter is highly efficient 
and is photon-noise-limited for point sources. Additional 
details about this polarimeter can be found in Magalhaes 
et al. imi . Pereyra IMiTii and Pereyra & Magalhaes 

imm . 

The data were collected on 20 and 21 June 1998. The 
measurements were made using a standard V filter with a 
1024x1024 CCD that covers an area of ~ 10’xlO’. Typical 
sequences of four, eight, or sixteen positions separated by 
22.5° of the half-waveplate were used depending on the 
magnitude of the star. Typical integration times for po¬ 
sition of the half-waveplate were a few seconds for the 
brighter stars and a few minutes for the dimmer stars. 
Our sample, along with a log of observations, is shown in 
Tabled The total sample is composed of sixteen objects 
and includes fourteen LRG and two NG. Th LRG were 
obtained from de la Reza & Drake (I1995H and Castillho 
et al. GS23- The object names are indicated in column 
(1) with their Li abundances in column (2). The galactic 
coordinates are in columns (3) and (4). The visual mag¬ 
nitudes are in column (5) and were obtained (when avail¬ 
able) from the literature. In particular, IRAS 19038-0026, 
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Fig. 3. Polarization maps for each target where foreground polarization 
was calculated. Each vector represents one object included in the average 
for each case (see Table |2j. The polarization scale is shown to the upper 
right of each map. The maps are centered on the coordinates of each 
target with its observed polarization also shown. The coordinates are 
2000 . 
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Table 2. Optical linear polarization measurements. 


Object 

Pobs 

®obs 

Pfor 

@ for 

N 

P/ap 

Pint 

®int 


(%) 

n 

(%) 

n 



(%) 

n 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

HD 90082 

0.272 (0.014) 

109.43 

0.649 (0.028) 

116.00 

9 

5 

0.389 (0.031) 

30.57 

HD 95799 

0.296 (0.021) 

96.03 

0.440 (0.014) 

111.93 

17 

5 

0.245 (0.025) 

41.72 

HD 96195 

1.056 (0.039) 

90.73 

1.127 (0.063) 

97.99 

10 

5 

0.285 (0.074) 

42.17 

HD 120602 

0.135 (0.030) 

70.23 







PDS 68 

0.305 (0.054) 

56.53 

0.104 (0.001) 

105.38 

32 

30 

0.335 (0.054) 

47.59 

HD 146850 

0.502 (0.041) 

102.13 







GCSS 557 

1.475 (0.066) 

43.03 

0.831 (0.001) 

58.37 

222 

20 

0.871 (0.066) 

28.46 

HD 176588 

0.663 (0.062) 

66.23 







IRAS 19012-0747 

1.903 (0.282) 

14.33 

1.786 (0.017) 

16.60 

26 

10 

0.187 (0.283) 

169.77 

IRAS 19038-0026 

0.418 (0.057) 

65.23 

0.351 (0.006) 

88.38 

97 

5 

0.309 (0.057) 

37.56 

HD 178168 

0.666 (0.022) 

59.43 

0.781 (0.018) 

47.56 

9 

5 

0.318 (0.028) 

108.87 

HD 112127 

0.032 (0.041) 

65.03 







V859 Aql 

2.074 (0.038) 

80.63 

0.989 (0.019) 

89.10 

38 

5 

1.164 (0.043) 

73.47 

HD 203251 

0.136 (0.054) 

143.13 







HD 190664 

0.655 (0.026) 

86.03 







HD 124649 

1.846 (0.021) 

63.23 

2.377 (0.076) 

68.99 

16 

5 

0.677 (0.079) 

175.48 


The polarization errors are indicated in parenthesis. 


with no previous V magnitude published, was identified 
as the S300232037086 object in the GSC 2.2 catalogue, 
and the photographic bands for this object were used in 
this compilation. The number of waveplate positions used 
in each object is indicated in column (6) and the integra¬ 
tion time used by waveplate position is in column (7). The 
radius of the aperture, which minimized the polarization 
error, used is shown in column (8). The date of each ob¬ 
servation is indicated in column (9) and identification for 
each object is shown in the last column. 

For the data reduction process, we follow the proce¬ 
dure indicated in Pereyra & Magalhaes (1200211 using the 
PCCDPACK package fPerevra l2()()()!> . This is a set of spe¬ 
cially developed IRAF scripts to study the polarization 
data in (eventually crowded) stellar fields. The linear ob¬ 
served polarization (P) and the polarization position angle 
(0, measured from north to east) for the sample are shown 
in columns (2) and (3) of Table [3 Corrections of the po¬ 
larization position angle to the equatorial system were ob¬ 
tained from polarized standard stars observed each night. 
Unpolarized standard stars were used to check the instru¬ 
mental polarization, which was found to be smaller than 
0.04%, so no correction for instrumental polarization was 
applied to the data. 

One star in our sample, LRG: HD 146850, is also 
present in the catalogue of stellar polarization by Heiles 
(MS- The Heiles values (P = 0.61±0.04%, 0 = 
89.1±1.6°) are consistent with our P value, but 6 presents 
a discrepancy that might indicate some variability. Our 
polarization measurements will be considered for the dis¬ 
cussion in the next section. 

Figured shows the histogram of observed polarization 
for the stars in our sample. Approximately half of the LRG 
stars have observed polarizations that are lower than 0.5%, 
but in a few cases high polarizations (>1.5%) are observed 
(IRAS 19012-0747 and V859 Aql), and one of the normal 


giants has an enhanced polarization value (~ 1.8%, HD 
124649). 

2.1. Foreground polarization 

The foreground polarization must be important if, for ex¬ 
ample, the line of sight is along the Galactic plane. Figure 
[2]shows the observed polarization for our sample as a func¬ 
tion of the Galactic latitude modulus. Clearly, for the ob¬ 
jects with \b | < 20°, the polarization covers a wide range. 
For higher latitudes, the observed polarization presents 
lower values (< 0.5%). Thus, for half of our sample the 
contribution of foreground polarization is important and 
must be subtracted. It is interesting to note that the ob¬ 
jects with the lowest observed polarizations (LRG: HD 
112127, HD 120602, and HD 203251) are at higher lat¬ 
itudes. We can conclude that the observed polarizations 
of these objects represent the very low or null foreground 
polarization in these directions. 

The foreground polarization for the objects in our sam¬ 
ple was estimated with the field stars in each CCD frame 
using PCCDPACK. This package provides the polariza¬ 
tion for each object present in a typical CCD frame, as 
well as the average Stokes parameters (Q and U) weighted 
by the errors for all the objects considered. With these pa¬ 
rameters the average foreground polarization can be cal¬ 
culated in a field of view that is very close to the target 
and just limited by the CCD size. We used a polarization 
signal-to-noise ratio ( P/crp ) of 5 or larger in order to se¬ 
lect objects to be included in the average (see Fig. HJ). We 
take this average as an estimate of the foreground polar¬ 
ization value. These are indicated in columns (4) and (5) 
of Table[2j along with the number of objects and the mini¬ 
mum threshold P/op used in each case in columns (6) and 
(7). For bright objects, which required short integration 
times, and/or for those with higher latitudes, it was usual 
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Table 3. Sample stellar parameters. 


Object 

rp a 

1 eff 

Log f 

[Fe/H]“ 

U-B b 

B-V 

V-R b 

R-P 

(B-V) 0 C 

E(B-V) 

Plx d 

Distance 

E(B - V) Bond 











(mas) 

(Pc) 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

(12) 

(13) 

HD 90082 

3600 

0.0 

0.0 

1.95 

1.67 

1.03 

1.15 

1.60 

0.07 

7.20 (4.80) 

139 r 5 ' 6 ° 

0.03 

HD 95799 

4900 

2.5 

0.0 

0.79 

1.00 

0.53 

0.51 

0.95 

0.05 

20.50 (11.50) 

49 —is 

0.01 

HD 96195 

3600 

-0.5 

0.0 

1.61 

2.28 

1.42 

1.56 

1.74 

0.54 

1.12 (0.74) 

893 + 1 3 7 5 3 5 9 

0.18 

HD 120602 

5000 

3.0 

-0.1 

0.59 

0.90 

0.68 

- 

0.90 

0.00 

8.09 (0.81) 

124 

0.00 

PDS 68 

4450 

2.6 

0.0 

1.90 

1.58 

0.88 

0.79 

1.14 

0.44 


- 


HD 146850 

4000 

1.5 

-0.3 

1.70 

1.48 

0.91 

0.79 

1.36 

0.12 

3.77 (0.85) 

265 /ll 

0.04 

GCSS 557 

3400 

0.0 

0.0 

- 

4.68 

3.13 

2.35 

1.57 

3.11 


3113 e 

0.56 

HD 176588 

4000 

1.5 

0.0 

2.10 

1.65 

0.97 

0.88 

1.38 

0.27 

4.25 (0.91) 

235 +*j 

0.05 

IRAS 19012-0747 

3800 

1.5 

0.0 

2.30 

1.83 

1.02 

0.93 

1.39 

0.44 


2330 e 

0.24 

IRAS 19038-0026 

3600 

1.0 

0.0 

- 

1.20 

- 

- 

1.35 

neg. 


- 


HD 178168 

4000 

1.0 

0.0 

2.17 

2.03 

1.09 

0.97 

1.43 

0.60 

17.70 (16.90) 

56 

0.01 

HD 112127 

4340 

2.1 

0.3 

1.43 

1.26 

0.90 

0.58 

1.26 

0.00 

8.09 (0.85) 

124 +« 

0.00 

V859 Aql 

4500 

2.5 

0.0 

1.44 

1.31 

0.82 

0.73 

1.12 

0.19 


- 


HD 203251 

4500 

3.0 

-0.3 

1.10 

1.22 

- 

- 

1.08 

0.14 

-0.20 (9.10) 

plx < 0 


HD 190664 

4650 

2.7 

0.0 

1.00 

1.16 

0.86 

- 

1.06 

0.10 

9.57 (0.96) 

104 +i2 
-10 

0.02 

HD 124649 

3750 

1.5 

0.0 

- 

1.70 

- 

- 

1.39 

0.31 

1.98 (0.97) 

505 

— 166 

0.09 


(from Castilho et a! 11999at . V859 Aql (from Reddy et al. 12002b HD 203251 (from Fekel & Balachandran 11053 i. and PDS68, HD 
190664, and HD 124649 (see text); (b) obtained from Castilho Ill999bt . except for: HD 120602, HD 112127, HD 203251, and HD 
190664 (from GCPC, Mermilliod et al. 1 1997 i. PDS 68 and V859 Aql (from Gregorio-Hetem et al. 119921 . HD 124649 from CDS. 


The B - V color index for GCSS 557 and IRAS 19038-0026 were obtained from the CSC 2.2 catalogue, using magnitudes in J and 
V photographic bands as B and V, respectively; (c) obtained from Bessell et al. f 19981 with interpolation (or extrapolation) of 
T e ff, log g and [Fe/H] (columns (2), (3), and (4)); (d) obtained from the Hipparcos catalogue, except for: HD 90082, HD 95799, 
HD 178168, and HD 203251 (from the Tycho catalogue). The parallax errors are indicated in parenthesis; (e) from Castilho et 

ai. cra¬ 


te find no field stars on the CCD frame. In these cases, 
it was not possible to obtain the foreground polarization, 
so those values are left blank in Table 0 In three of these 
cases (LRG: HD 112127, HD 120602, and HD 203252), the 
observed polarization has a P/ap lower than 5 and these 
objects can be considered unpolarized. These objects are 
those located at higher latitudes as mentioned above. 


2.2. Intrinsic polarization 

For those objects with an estimated foreground polariza¬ 
tion (Table0, it was possible to obtain the intrinsic polar¬ 
ization. As the foreground polarization is an additive com¬ 
ponent included in the observed polarization, the intrinsic 
Stokes parameters (Q = Pcos(20) and U = Psin(20)) are 
as follows: 


Qint — Qobs ” Qf or 
Pint — Pobs ” Ufor- 


The intrinsic polarization (Pint) and its polarization 
angle (0i n t) are obtained from: 

Pint = (Qint + Pint 
®int 2 aiCtan( Pint / Qint ) ■ 



The intrinsic polarization estimated for ten objects 
(9 LRG and 1 NG) are indicated in columns (8) and 
(9) of Table 0 and shown in Fig. 0j. In general, the 
intrinsic polarization levels are significant. Seven LRG 
(five of them with P/ap > 5) have intrinsic polariza¬ 
tions between (0.19-0.34)%, but in some cases high values 
(> 0.5%) are obtained (LRG: GCSS 557 and V859 Aql; 
and NG: HD124649). As an extreme case, one of the ob¬ 
jects with high observed polarization (IRAS 19012-0747, 
1.90%) presents negligible intrinsic polarization within 5er 
(0.19%). 

A quantitative analysis of intrinsic polarization is be¬ 
yond the scope of this paper; nevertheless, the polariza¬ 
tion models of circumstellar dust shells from Johnson and 
Jones UTiMl) are useful for investigating the geometry of 
the shells in our sample. In that work, the shells around 
evolved stars were modeled as ellipsoids of revolution and 
the results pointed out the importance of taking into ac¬ 
count the inclination of the ellipsoid to the plane of the 
sky (i’). For an ellipsoid of any given axial ratio, the incli¬ 
nation would decrease the polarization by ~ cos 2 (P) (or 
~ sin 2 (?), where i is the inclination to the line of sight). 
Therefore, a tilted ellipsoid would need a higher axial ratio 
to produce the same amount of polarization as a nontilted 
ellipsoid would. Comparing the red giant edge-on mod¬ 
els in Johnson and Jones |MD, polarization of ~ 0.3% 
was found for slightly prolate optically thin shells (as in 
R Vir) and ~ 2.3% for prolate optically thick shells (as 
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in RU Vir). This last case is more consistent with the full 
range of intrinsic polarization in our sample [0.19—1.16]%. 
Assuming the dependence with sin 2 (z), our sample is best 
represented by pole-on objects (i < 30°), except for GCSS 
577 and V859 Aql. These two objects have a high chance 
of being closer to edge-on (probably, i > 45°). A more de¬ 
tailed modeling using Monte Carlo methods developed by 
our group (Carciofi et al. 120041) is planned. 

3. Discussion 

3.1. E(B - V) and distances 

Correlations between the observed polarizations and the 
color excess can also help to investigate if an intrinsic color 
excess is associated with intrinsic polarization. This can 
be especially important for objects with undetermined in¬ 
trinsic polarization. For this purpose, we computed the 
color excess, E(B - V), for the objects in our sample, using 
the information in the literature when available (see Table 
0 ). 

The effective temperatures, surface gravities, and 
metallicities are shown in columns (2), (3), and (4) of 
Table 0 Three objects (PDS 68, HD 190664, and HD 
124649) have unknown stellar parameters to date and 
the listed values were calculated following Castilho et al. 
(|200Qll using the near infrared colors (see Sect. 13.2. Tl) . 

The color indexes U -B, B -V, V -R, and R -I, when 
available, are shown in columns (5) to (9) of TableEH along 
with the intrinsic color index ( B - V)q. In particular, GCSS 
557, with no previously published B - V index, was iden¬ 
tified as S300112223323 in the GSC 2.2 catalogue and its 
B - V color index compiled here was constructed using the 
available photographic bands ( J and V, respectively). As a 
comparison, the visual magnitude V=13.26 from Castilho 
( 1999b ) agrees with the V=13.50 photographic band mag¬ 
nitude (from GSC 2.2). An important discrepancy is found 
between its color ( V-R = 3.13) from Castilho (I1999bl) and 
the (V -F = 0.64) index using the red photographic band 
( F ) from GSC 2.2. Then, the B - V index used here for 
this object must be taken with care. Finally, the calcu¬ 
lated color excess, E(B - V), for our sample is indicated in 
column (10) of Table 0 Negative values were not consid¬ 
ered, as in case of IRAS 19038-0026. 

Figure 2k shows the correlation between the observed 
polarization and the color excess, E(B -V). The upper 
limit for optimum alignment of polarizing dust grains in 
the diffuse ISM is also shown (Py / E{B - V) = 9% mag -1 , 
Serkowski et al. mraa . All the objects, except V859 Aql 
and HD 120602, are under this limit, so we might con¬ 
clude that the observed polarization has an interstellar 
origin. Nevertheless, the non-negligible polarization found 
in Sect. l2~2l tell us that a fraction of the excess color must 
originate in the circumstellar material in these stars. HD 
120602 has an E{B - V) = 0, and its position in Fig. 2k 
probably reflects the very low or null foreground polariza¬ 
tion in this direction as we noted in Sect. im In contrast, 
the case of V859 Aql could indicate a very important con- 
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Fig. 4. Correlation between observed polarizations and 
color excess (a) using intrinsic colors, and (b) using the 
Bond extinction law. The black dots are the LGR, and the 
white dots are the NG. The dashed line represents the up¬ 
per limit for optimum alignment for polarizing dust grains 
in the diffuse ISM (Py /E(B -V) = 9% mag -1 , Serkowski 
et al. ITO7BI . The box inside (b) figure is a zoom to show 
more detail. The letters indicate the ID for each object as 
it appears in column 10 of Table |T] 


tribution of the circumstellar envelope to the excess color 
that is observed. Both V859 Aql and GCSS 557 are the 
two LRG with the highest computed intrinsic polarization, 
and they have a high chance of presenting an asymmetric 
spatial distribution of circumstellar dust (envelope) prob¬ 
ably viewed edge-on. 

An alternative method to obtain the color excess is 
to apply the Bond ITMH extinction law, but a distance 
estimate is needed. Distances were derived from paral¬ 
laxes given in the Hipparcos and Tycho catalogues when 
available. The parallaxes and distances are indicated in 
columns (11) and (12) of Table 21 respectively, but neg¬ 
atives values were not considered (HD 203251). Objects 
with parallax errors that produce distance errors that are 
higher than 100 pc must be considered with care (HD 
90082, HD 96195, HD 178168, and HD 124649). For GCSS 




















Pereyra et al.: Polarimety of Li-rich giants 


Table 4. Infrared data. 


Object 

2MASS ID 

J 

H 

K 

Flag“ 

(J-H) o 

6 (H-K) o 6 

E(J-H) 

E(H-K) 

IRAS ID 



(mag) 

(mag) 

(mag) 







(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

HD 90082 

19062377-0021222 

5.43 

4.48 

3.98 

ADD 

0.94 

0.20 

0.02 

0.30 

10204-6135 

HD 95799 

11022509-5845373 

6.21 

5.71 

5.55 

AAA 

0.54 

0.04 

-0.04 

0.12 

none 

HD 96195 

11042916-6257439 

3.48 

2.43 

2.07 

DDD 

0.91 

0.22 

0.15 

0.14 

11024-6241 

HD 120602 

13502469+0529498 

4.74 

4.22 

4.10 

DDD 

0.51 

0.03 

0.01 

0.09 

13479+0544 

PDS 68 

13565727-4208098 

9.90 

9.16 

8.89 

AAA 

0.68 

0.07 

0.06 

0.20 

13539-4153 

HD 146850 

16190043-1452220 

3.27 

2.46 

2.27 

DDD 

0.87 

0.11 

-0.06 

0.08 

16161-1445 

GCSS 557 

18265916-1441487 

4.86 

3.34 

2.71 

ECD 

1.03 

0.23 

0.50 

0.40 

18241-1443 

HD 176588 

19010951-0426201 

4.10 

3.19 

2.96 

DDD 

0.84 

0.11 

0.07 

0.12 

18585-0430 

IRAS 19012-0747 

19035796-0742345 

6.73 

5.63 

5.16 

AAA 

0.92 

0.14 

0.18 

0.33 

19012-0747 

IRAS 19038-0026 

19062377-0021222 

5.43 

4.48 

3.98 

ADD 

0.98 

0.17 

-0.02 

0.33 

19038-0026 

HD 178168 

10220747-6151005 

4.00 

3.03 

2.79 

DDD 

0.82 

0.12 

0.15 

0.12 

19049-0234 

HD 112127 

12535573+2646479 

5.00 

4.34 

4.12 

DCE 

0.70 

0.09 

-0.04 

0.14 

12514+2703 

V859 Aql 

19310123+0523533 

7.57 

6.82 

6.60 

AAA 

0.66 

0.06 

0.09 

0.15 

19285+0517 

HD 203251 

21212967-1509216 

5.83 

5.18 

5.02 

AAA 

0.67 

0.06 

-0.02 

0.10 

21187-1522 

HD 190664 

20061222-0404414 

4.51 

3.95 

3.85 

ECC 

0.62 

0.06 

-0.06 

0.03 

20035-0413 

HD 124649 

14162210-5319365 

4.84 

3.78 

3.65 

DDD 

0.92 

0.14 

0.14 

-0.01 

14130-5305 


Notes: (a) 2MASS JHK photometric quality flag (from ”A”: very good, to ”E”: very poor); (b) from Bessell et al. nr 



Fig. 5. Correlation between observed polarizations and 
distance. The black dots are the LGR, and the white dots 
are the NG. The box inside the figure is a zoom to show 
more detail. The letters indicate the ID for each object as 
it appears in column 10 of Table |T| 


557 and IRAS 19012-0747, with unknown parallaxes, the 
distances were obtained from Castilho et al. < 20001 , who 
use a color-magnitude diagram for Hipparcos field stars 
(Perryman et al. 1199511 . We must note here that three 
objects (HD 90082, HD 95799, and HD 178168) with dis¬ 
tances from Tycho catalogue show significant discrepan¬ 
cies with the (apparently overestimated) values found by 
Castilho et al. 1200011 for these stars. In the following, our 
tabulated distances to these 3 stars will be used. 

Figure0plots the observed polarization and distance. 
Seven LRG and one NG are located within the first 300 
pc with observed polarizations lower than 0.75%. It is in¬ 
teresting to note that GCSS 557 with its high observed 
polarization (and also intrinsic polarizations) is located at 


an extreme distance. This is consistent with the high level 
of foreground polarization calculated in this direction. The 
important foreground contribution in the line of sight to 
IRAS 19012-0747 (see Sect. 12.2|) is also consistent with the 
larger assumed distance (2330 pc) for this object and with 
the fact that it is located at low galactic latitude. 

With distance information we can obtain the color ex¬ 
cess using the Bond dMj extinction law, E(B -V)Bond- 
For objects with b > 60°, E(B -V)Bond = 0; for b < - 
60°, E(B -V) B ond = 0.03; and for |6| < 60°, E(B -V) Bon d 
= 0.03csc6[l-exp(-0.008rsin6)] where r is the distance in 
parsecs and b is the Galactic latitude (in absolute value). 
This is shown in column (13) of Tabled For GCSS 557 and 
IRAS 19012-0747, the results are consistent with Castilho 
et al. m . If we consider that E(B -V) includes the 
contribution of the foreground ISM along with an intrin¬ 
sic color excess (such as in Fig. Uk), and if E{B - V) Bo nd 
represents just the contribution of foreground ISM to a 
given line of sight, we would expect that an intrinsic color 
excess could explain the intrinsic polarization observed. 
Fig. EH> shows the correlation between the observed polar¬ 
ization and the color excess using the Bond law. We see 
that a large part of the LRG present an excess of polariza¬ 
tion when we compare with the upper limit of polarization 
with an ISM origin (Serkowski et al.: 1 9751 ). 

Just three LRG with non-negligible intrinsic polariza¬ 
tions (HD 96195, GCSS 557, and IRAS 19012-0747) are 
under the Serkowski limit, but in two cases (GCSS 557 and 
IRAS 19012-0747) the inferred distance may have been 
overestimated as mentioned earlier. If it is the case, E(B 
- V)Bond for these objects probably has a lower level and 
therefore these points are located more to the left. The 
distance error for HD96195 obtained from its parallax er¬ 
ror is bigger than 100 pc and also must be taken with 
care. 
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Fig. 6. The (J - H) x (H - K) color-color diagram. The 
locus of the giant branch is shown by the solid line with 
spectral types indicated from Koornneef lfT3%31) and the 
dashed lines follow the reddening vectors taken from Rieke 
& Lebofsky Ill985|l . The black dots are the LRG, and the 
white dots are the NG. The letters indicate the ID for each 
object as it appears in column 10 of Tabled 


Seven LRG and two NG have an observed polarization 
above the Serkowski limit, and in four of them (LRG: HD 
90082, HD 95799, and HD 178168; and NG: HD 124649) 
we detected non-negligible intrinsic polarization. The two 
LRG with the lowest observed polarizations (HD 120602 
and HD 112127) are also slightly above this limit, but 
in both cases E(B -V)Bond = 0 because |b| > 60°. This 
is consistent with our assumption of unpolarized objects 
given in Sect.Q 

It is interesting to note that three objects (LRG: HD 
146850, HD 176588; NG: HD 190604) with significant ob¬ 
served polarization and without estimated foreground po¬ 
larization also appear above the Serkowski limit in Fig. 
Q]). Thus, the correlation P 0 b s vs. E(B -Vjsond also gives 
information on a possible intrinsic color excess when an 
estimate of intrinsic polarization is not possible. 

In general, analysis of the correlations of P 0 bs with 
the Bond law supports the idea that the non-negligible 
intrinsic polarization observed in LRG (and also in NG) 
has an origin in the circumstellar material. 

3.2. Correlations with infrared colors 

In a simple scenario, circumstellar dust in LRG must ab¬ 
sorb light from the central object and radiate thermally 
in the infrared. Then, dust scattering will produce po¬ 
larization, and a correlation between the optical intrinsic 
polarization and the excess (near or far) infrared emis¬ 
sion is expected. Additionally, these correlations can give 
us information about the type of scatterer in the non- 
spherically symmetric envelope. 
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Fig. 9. (a) Intrinsic polarization versus [25 - 12] IRAS 
color; (b) Intrinsic polarization versus [60 - 25] IRAS color. 
The black dots are the LRG, and the white dots are the 
NG. The letters indicate the ID for each object as it ap¬ 
pears in column 10 of Tabled 

3.2.1. Near infrared colors 

We used the 2MASS All-Sky Catalog of Point Sources 
(Cutri et al. 1 2003 1 to identify the NIR colors for the ob¬ 
jects in our sample. The results are indicated in Tabled 
In columns (1) and (2), we show the named object and its 
NIR counterpart identified by the 2MASS ID, respectively. 
The J, H 1 and K colors are shown in columns (3), (4), and 
(5). Column (6) indicates the 2MASS JHK photometric 
quality flag for each measurement (from ”A”: very good, 
to ”E”: very poor). 

In Fig. El we present the J - K versus H - K diagram 
for our sample, together with the positions of the giant 
branch IKoornneef 11983 ) and the directions of the red¬ 
dening vectors (Rieke & Lebofsky 11985 1 between G8III 
and M4III types. The majority of the objects have the 
colors of reddened photospheres, but two of LRG stars 
(HD 90082 and IRAS 19038-0026) appear to the right of 
the reddening line for G8 giants, showing a slight excess 
at 2.2 /jm. With the NIR intrinsic colors obtained from 
Bessell et al. (119981 (cols. (7) and (8) in Table EJ, we cal¬ 


culated the NIR excess colors E(J - H) and E(H- K) (cols. 
(9) and (10) in Table EJ for our sample using the stellar 
parameters from columns (2), (3), and (4) of Table El 

The correlations between the intrinsic polarization and 
the NIR excess color are shown in Fig. [3 In order to 
serve as an objective criterion to quantify the correlation 
between Pi nt and the NIR excess color, we choose the 
Spearman rank-order correlation coefficient (r s ), to which 
a statistical significance can be attached. The calculated 
value of r s for the Pi nt versus E(J - H) relation (Fig. EH) 
is 0.10 with a 78% probability, where a low probability 
means high significance. For the Pint versus E(H - K) re¬ 
lation (Fig. Q d), r s is 0.07 with an 85% probability. With 
this information, we can conclude that there is no statisti¬ 
cally significant correlation between and NIR excess 
color. V859 Aql and GCSS 557 with the highest intrinsic 
polarizations present low and high NIR excess color val¬ 
ues, respectively; while HD 90082 and IRAS 19038-0026 
with an excess in 2.2 /iin (Fig. El have lower intrinsic po¬ 
larization values. Using the NIR results, we can conclude 
that NIR emission either (a) is not responsible for the 
optical Pint in LRG or (b) is responsible for the optical 
Pint in LRG and the majority of LRG in our sample are 
pole-on objects except V859 Aql and GCSS 557. 

3.2.2. IRAS colors 

In order to investigate if the large far infrared excess for 
LRG stars observed in IRAS colors is correlated with the 
optical intrinsic polarization, we identified each object of 
our sample with its IRAS counterpart (IRAS Point Source 
Catalogue). The results are shown in column (11) of Table 
El and just one object (LRG: HD 95779) of our sample 
does not present an IRAS counterpart. 

Figure El shows the distribution of our sample in the 
IRAS 12-25-60 (im color-color plane. The normal giants 
HD190664 and HD 124649 are coincident with the lower 
side of LRG stars. The locus of LRG stars in our sample is 
consistent with groups 1 and 3 from Gregorio-Hetem et al. 
1199811 . FigurcElplots the optical Pint versus [25 - 12] and 
[60 - 25] IRAS colors. Considering the objects with higher 
intrinsic polarizations (GCSS 557 and V859 Aql), an ap¬ 
parent correlation is observed with [25-12] color and an 
anticorrelation with [60 -25] color. It could be consistent 
with our assumption of a favorable (edge-on) viewing an¬ 
gle for the envelope in these two objects. Again, the lower 
polarization observed in the rest of LRG stars would be 
reflecting objects with an intermediary viewing angle or 
pole-on envelopes. To quantity this finding, we calculated 
the r s coefficient (as in Sect. rrm for all the data in Figs. 
EH^b- The calculated value of r s for the Pi nt versus [25 - 
12] color relation is 0.15 with a 70% probability; and for 
the Pint versus [60-25] color relation, r s is -0.72 with a 3% 
probability. The positive (and less significant) and nega¬ 
tive (and more significant) correlations with IRAS colors 
also can imply that dust with preferential emission in ~ 
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Fig. 10. Spectral energy distributions. In squares, the optical data; in triangles, the NIR 2MASS data; in diamonds, 
the MSX data (band quality > 0); and in circles, tha IRAS data. In solid lines, the two black-body fit to the data. 
In dashed lines, the black-body with T e ff ; and in dotted lines, the black body with the calculated dust temperature 
( Tdust )• 


25 ;tm is responsible for the optical intrinsic polarization 
observed. 

3.3. SEDs 

In order to check the typical dust temperature (Tdust) for 
the envelopes that are responsible for the optical polariza¬ 
tion, we proceeded to construct the spectral energy distri¬ 
bution (SED) for the objects in our sample. For that we 
used the optical, near infrared, and far infrared data avail¬ 


able (Tables Cl El and n. We also used, when available, 
the intermediary infrared data from the MSX6C Infrared 
Point Source Catalog (Egan et al. (20031 ). The objects with 
MSX data are indicated in column (2) of Tabled The flux 
calibration was made using Cox cnnni for the optical data 
and Cohen et al. 1200811 for the NIR data. 

The results are shown in Fig.^3 The optical and NIR 
data were not corrected by reddening. We also show an ad¬ 
ditive two black-body fit that includes a component with 
a fixed effective temperature (column (2) in Tabic |3j) rep- 
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Table 5. Additional information. 


Object 

(1) 

MSX6C ID 

(2) 

E-h s ' 

HD 90082 

G286.2685-03.9091 

338 

HD 95799 

G289.1898+01.1718 

- 

HD 96195 

G291.1180-02.5694 

210 

HD 120602 

- 

- 

PDS 68 

- 

191 

HD 146850 

- 

239 

GCSS 557 

G017.0984-01.3757 

240 

HD 176588 

G030.0809-04.1814 

205 

IRAS 19012-0747 

- 

215 

IRAS 19038-0026 

G034.3214-03.4902 

257 

HD 178168 

G032.5404-04.7223 

105 

HD 112127 

- 

609 

V859 Aql 

- 

217 

HD 203251 

- 

- 

HD 190664 

- 

221 

HD 124649 

- 

246 


resenting the central object and a second component with 
the temperature as a free parameter representing Tdust 
for the envelope. The fit was done over the NIR, MSX, 
and IRAS data to minimize the effect of reddening in the 
optical bands. The IRAS higher limits were not consid¬ 
ered in the fit, while three objects (HD 96195, HD 146850, 
and HD 124649) with apparently good quality IRAS Tjoo 
present cirrus contamination {cirr3 /Fqq > (1- 5), Ivezic & 
Elitzur fl995ll . and Fioo also were not considered. Typical 
dust temperatures for the envelopes were between 190 K 
and 260 K (see column (3) in Table 0 associated with 
maximum emission between 19 //m and 14 /mi, respec¬ 
tively. This fact agrees with our conclusion about grains 
emitting in ~ 25 /mi as responsible for the intrinsic polar¬ 
izations (see Sect. I3T5I . As a comparison, Reddy et al. 
( 2002 ) find T dust ~ 250 K for V859 Aql, which compares 
with our value ( T du at = 217 K) very well. 

4. Conclusions 

Optical polarimetry for a sample of fourteen LRG stars 
was obtained. For nine of them intrinsic polarization was 
estimated using field stars. Seven LRG (five of them with 
P/ap > 5) have lower but non-negligible intrinsic polar¬ 
ization (0.19 — 0.34)%, and in two cases (V859 Aql and 
GCSS 557) intrinsic polarizations higher than 0.5% are 
found. These results indicate that an asymmetric spatial 
distribution of circumstellar dust is present in LRG (but 
probably also in normal giants). An excess in observed 
polarization when it is correlated with the optical excess 
color gives additional support to the circumstellar origin 
of the intrinsic polarization in LRG. The optical intrin¬ 
sic polarization in LRG is not correlated with the near 
infrared excess but is correlated with far infrared emis¬ 
sion. This would suggest that grains emitting in ~ 25 /tm 
are responsible for the optical intrinsic polarization, and 
the higher intrinsic polarization levels would indicate a 


favorable (edge-on) viewing angle for the envelopes, as in 
V859 Aql and GCSS 557. Analysis of spectral energy dis¬ 
tributions for the sample provides an estimate of the dust 
temperature for the envelopes, which are mainly between 
190 K and 260 K. Our findings indicate that non-spherical 
symmetries may appear as early as the RG phases of stel¬ 
lar evolution. 
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